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We analyze the spherically symmetric Einstein field equation with a massless complex scalar field. 
We can use the Newtonian solutions to fit the rotation curve data of spiral and dwarf galaxies. From 
the general relativistic solutions, we can derive high gravitational redshift values. 
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I. INTRODUCTION 

The rotation curves for galaxies or galaxy clusters 
should show a Keplerian decrease v ~ \/V^ a ^ ^ ne P om t 
where the luminous matter ends. Instead one observes 
flat rotation curves beyond the galaxies jjj . A linear ra- 
dial increase of the mass function of galaxies and galaxy 
clusters have been derived from these observations [g): 
x. Several models have been discussed where 



limit 



M 

either non-Newtonian gravity |3J or non-interacting mat- 
ter, dark matter Q, are introduced to solve this prob- 
lem. For several classes of gravitational theories, it was 
recently shown that the introduction of dark matter is 
necessary Hj; cf. ||. Massive compact halo objects, so- 
called MACHOs, consisting of baryonic matter are also 
not able to solve this problem Qj. 

In 1933 Zwicky || was the first who suggested the ex- 
istence of dark matter in galaxy clusters by investigating 
the Coma cluster. The total mass needed to gravitation- 
ally bind this cluster exceeds the amount of the luminous 
matter by roughly an order of magnitude. Three years 
later Smith proved this for the Virgo cluster || . Begin- 
ning of the 1970's one was able to extend the measure- 
ments of the rotation curves of galaxies so that higher 
mass to luminosity relations could be found: after some 
radius the rotation curves revealed that there is more 
mass than contained within the luminous matter [ p~0| . 
The explanation of a linearly increasing mass was first 
given by Freeman ]TT[ | providing a spherical halo. The 
investigations for determining the radius of a dark mat- 
ter halo have to go beyond the HI measurements of the 
1980's fil| , e.g. by studying satellite galaxies Jli|, using 
the weak lensing of background galaxies by foreground 
dark haloslflj], or looking into quasar absorption lines 
p5| ; cf. ]16[ . From these investigations, halo radii of 
more than 200kpc are inferred, for our Galaxy 230kpc 
p7[ , and recent results from satellite galaxies of a set of 
spiral galaxies show even more than 400kpc |18fl . Re- 
cently, measurements of rotation curves of high redshift 
galaxies have been carried out [jl9| . 

We present a solution class of a massless complex scalar 
field minimally coupled to the Einstein equation p0[ . For 



the Newtonian types of these solutions, we can fit rota- 
tion curve data of spiral and dwarf galaxies. The limiting 
value of the orbital velocity is determined by the central 
amplitude of the scalar field. The frequency of the scalar 
field determines the halo characteristics of mass and den- 
sity near the center. 

For the general-relativistic (GR) solutions, we can 
show that they provide large gravitational redshifts. We 
discuss how emission and absorption lines produced in 
the highly relativistic potential of our model can be un- 
derstood. Near the center of these solutions, we find 
rotation velocities of about 10 5 km/s so that high lumi- 
nosities can be expected. Spacetime singularities do not 
appear within these GR solutions. 

A self-gravitationally massive com plex scalar field is 
utilized for the so-called boson stars p4|| . These bo- 
son stars have no physical singularities as in our solutions 
or in the case of neutron stars. However, real massless 
p5[ or massive [^6| scalar fields (in each case without con- 
served Noether current) cannot prevent the formation of 
a singularity (cf. also the exact solution in ^7j). Such 
behavior of singular solutions is supported by analytical 
investigations of Christodoulou |28| and numerical calcu- 
lations of Choptuik p9] . 



II. EINSTEIN-SCALAR-FIELD EQUATIONS 

The Lagrange density of a massless complex self- 
gravitating scalar field reads 



C = 



1 



1 



-R + g^(d^*)(d^) 



(1) 



where R is the curvature scalar, k = 8nG, G the grav- 
itation constant (h = c = 1), g the determinant of the 
metric g^u , and $ the massless complex scalar field. Then 
we find the coupled system 



-kT m „($) 



□ $ = , 



(2) 
(3) 



where 



1 



V = (d^*)(d u $) - -9^[g™(d a $*)(d K $)} (4) 



is the energy-momentum tensor and 



a = d„ 



VTWdv /V\g\ 



(5) 



the generally covariant d'Alembertian. 

For spherically symmetric solutions we use the follow- 
ing static line element 



ds 2 = e v ^dt 2 - e x ^dr 2 - r 2 {M 2 
and, for the scalar field, the ansatz 
$(r,i) = P(r)e- iwt , 



+ sin 2 i? dip 2 ) (6) 



(7) 



where w is the frequency of the scalar field. 

The non-vanishing components of the energy-momen- 
tum tensor are 



T Q =p=-T 1 1 = Pr 

= X -[u?P 2 {r)e- v + P' 2 (r)e- X } 
T-2 2 = T 3 3 = -p ± 



1, 



(8) 
(9) 



where ' = d/dr. As equation of state, we find p = p r — 
p ± + P' 2 {r)e- X . 

The decisive non- vanishing components of the Einstein 
equation are 



v' + A' = n(p + Pr)ve > 



A' 



Kpre 



-(e A ~l) 



(10) 
(11) 



and two further identical components which are fulfilled 
because of the Bianchi identities. 

The differential equation for the scalar field is 



P"{r) + 



A' 



- ) /"(/■') +, x ■'\'-P[r) = . ( 1.2) 



A typical behavior of the scalar field of Newtonian kind 
is demonstrated in Fig. |l]; the metric potentials are al- 
most constant so that we do not show them here; cf. Sec- 
tion [II. General relativistic solutions are shown in Sec- 
tion [X. 

For the rest of our paper, we employ the redefined 
quantities x :— Lur and a := \J k/2 P. The numerical cal- 
culation was carried out by using a Runge-Kutta-routine 
of the Fortran- libraries IMSL/NAg. In order to get reg- 
ular solutions at the origin for the system of differential 
equations jlO|), (JTT|), and (|l2|), we have to put the initial 
conditions cr'(0) = and A(0) = 0. 

The system (|l0|), (|ll|), and jl2| ) possesses the self- 
similarity: x — > kx, A — > A, e v — > fc 2 e I/ , and a — > a. 
This self-similarity means that one has only one free pa- 
rameter, namely the initial value of the scalar field. One 
scales the solution simply with the second initial value 
for v. 
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FIG. 1. The massless complex scalar field a = y k/2 P 
depending on the dimensionless coordinate x = cor with the 
initial values a (0) = 3.8xl(T 6 [a(0)(x op t)] and u(0) = 0. This 
solution belongs to the rotation curve with Mi = —18.29 in 

Fig. a 



III. NEWTONIAN SOLUTION 

The Newtonian solutions of our Lagrangian (l|) are 
characterised by almost constant metric potentials v, A. 
The scalar field equation ( |l2] ) can then be rewritten into 
the following form: 



2a' 



+ (7 = 



(' = d/dx) which has the solution 

u{x) = — [Asin(x) + Bcos(x) 

x 



(13) 



(14) 



where A, B are some constants. For B ^ 0, the solution 
is singular at the origin, why we can rule out it. 
The Newtonian form of the energy density reads 



A 2 

P(x) = — 



1 - 



sin(2a;) sin 2 (a;) 



(15) 



For small x, we have p oc A 2 [l — 2x 2 /9 + x 4 /45], i.e. it 
is proportional to a constant. For dwa rf g alaxies, one 
observes such a behavior of the density [f50[ where good 
fitting results are obtained by using the empirical isother- 
mal density profile p(x) oc l/(x 2 + x 2 ), where x c is the 
core radius. Hence, our constant A resembles a core ra- 
dius. By applying a maximum halo model, i.e. with- 
out baryonic matter, the value of A is higher than for 
a model including baryonic matter. Therefore the core 
radius increases if one adds baryonic matter as found in 

a cf - & 

The general solution of Eq. ( pT| ) is 
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FIG. 2. The mass function in units of [a 2 (0)(i op t)/(a«t)] 
increases near the center like x 3 and at higher radial values 
linearly. We use here the same initial values as in Fig. [jj so 
that this curve shows the behavior of mass for the rotation 
curve with Mi = -18.29 in Fig. §. 



X 

with the mass function M(x) — 
the Newtonian formula (cf. Fig. 



(16) 

fop(Ot 2 d(. We find 



M(x) = A 2 { x + 



cos(2x) — 1 
2x 



(17) 



Comparing with the expected Newtonian result M — 
vf imit x, we see that vu m it — A; hence, the amplitude 
of the scalar field at the origin determines the limiting 
orbital velocity. For small x, the mass function behaves 
like 



M(x) oc A 2 x a /3 - 2A z x°/A5 + 0(x') 



(18) 



This corresponds to the constant density at the center. 
At higher radial distances, the mass function shows the 
linear behavior. 

Asymptotically, following (|l6|), the metric potential 
e~ A approaches the value C 2 := 1 — A 2 , where C 2 < 1. 
After a redefinition of the coordinate x — > C~ Y x , the 
asymptotic space has a deficit solid angle. The area of 
a sphere of radius x is not Airx 2 , but 4irC 2 x 2 ; cf., e.g., 
analogous results for global monopoles and global tex- 
tures |Q where one also finds a lin ear i ncre ase o f the 
mass function. We show in Sections VII and VIII how 



one can find a closure of the solutions and avoids this 
asymptotic problem. 

Following (|l0|), the second metric potential behaves 
asymptotically like e v — > x K , where K = 2A 2 /C 2 = 
2(1/C 2 — 1) > 0. The behavior of both metric potentials 
could be confirmed numerically. 

A question is the validity of the formulas (fl4|), (15), 
and @. We found for cr(0) = 1(T 3 and v(0) = at 



x = 10 a deviation of 0.001%, at x = 100 of 0.3% and 
at x = 1000 of about 4% . This shows clearly that these 
formulas should be used only near the center. All figures 
in this paper were produced by using the numerical so- 
lutions. If one is interested into the asymptotic behavior 
of the solutions, the formulas can still be used asymptot- 
ically (e.g. for the calculation in the Tables) because one 
can still derive the order of magnitude from them as we 
confirmed numerically. For higher intial values of a, only 
numerically determined solutions can be used. 



IV. ROTATION CURVES 

In this Section we shall model rotation curves of dwarf 
and spiral galaxies. Observations show that rotation 
curves are becoming flat (constant orbital velocity) in the 
surrounding region of galaxies where data are received 
from the 21cm wavelength of neutral hydrogen (HI). But 
the mass density of the neutral hydrogen is not sufficient 
to explain this velocity behavior. Therefore, we introduce 
dark matter consisting of a massless complex scalar field 
which interacts with the luminous matter exclusively by 
the gravitational force. We apply only Newtonian solu- 
tions of our model. 

For the static spherically symmetric metric (^) consid- 
ered here, circular orbit geodesies obey 



1 



1 



2 

M(r) 



\n P rr 2 e x+ » 



~n Pr r 2 e X+ » 



(19) 



which reduces outside of matter for a weakly gravita- 
tional field into the Newtonian form i> 2 Newt = M (r) jr. 
But, for p r ^ 0, we have to use the general-relativistic 
formula. By using the Newtonian solutions of Section III 
within (||), we find 



K = A 2 



1 - 



sin(2a;) 
2x 



cf. a generic rotation curve in Fig. |^. 
Asymptotically, from (|l9|), we have 



■A 2 / 2] 



v 2 v = e v [A 2 /2- 



(20) 



(21) 



what means that both the first Newtonian part e v (e x — 
l)/2 and the second matter part contribute the same 
amount to the rotation velocity. For e v = 1, we have 
v v = A as we had already derived from the general New- 
tonian solution for the mass function. 

The next step is now to take observational data for spi- 
ral and dwarf galaxies and try to model them by using our 
model and a model which describes the luminous matter 
distribution. For spiral galaxies, we use both the univer- 
sal rotation curves of Persic, Salucci, and Stel |3j| and 
some individual ones. Persic et al. confirmed by investi- 
gating data of 967 spirals that the structural properties of 
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FIG. 3. A generic rotation curve for an initial value of the 
scalar field cr(0) = A = 3 x 10" 3 . The Newtonian (M/x) and 
the pressure (p x ) contributions are shown separately. In the 
rotation curve fits, mainly the first part up to the maximum 
is used. The rotation curve oscillates around the asymptotic 
value of 90km/s. The velocity v is measured in [km/s] while 
the dimension of the radial coordinate r depends on the choice 
of ui\ x is dimensionless. 



dark and visible matter are linked together. This means 
that a spiral galaxy with low luminosity is stronger dom- 
inated by a dark matter halo than a spiral one with high 
luminosity. This is also one possible statement of the 
Tully-Fisher relation Q]. From this, it follows that a low 
luminosity spiral has a rather increasing rotation curve 
and a high luminosity spiral galaxy a rather decreasing 
rotation curve. 

In the following, we model the universal rotation curves 
by a combination of a stellar disk and our halo. The rota- 
tion curve for the stellar disk follows from an exponential 
thin disk light distribution fl33|| 



Vdisk( x )/ V ( X opt) 

0.325 



0.72 



2.5 



(M* B - M B ) 



1.97a: 1 - 22 



0.78 2 ) 1 - 43 



(22) 



where M B = 



20.5 is the absolute magnitude in the 
blue band (corresponding to logi* = 10.4) and Mb = 
-0.38 + 0.92Af/ (Mi from the I band); this formula can 
be used within the range 0.04 ~ x/x op t < 2. The x 
depending factor arrives from approximations of modified 
Bessel functions (see below); the constants in the Mb 
depending factor are arranged that it gives the best fit for 
our universal rotation curves; cf. j33|. The total rotation 
curves result then from 



Vtotal 



disk 



halo 



(23) 



The outcome can be seen in Fig. |j. One recognizes a very 
good agreement of the data. 

From this fit, we find that the amount of disk matter 
has to decrease with decreasing luminosity, i.e. decreasing 
absolute magnitude. This is a consequence of the Tully- 
Fisher relation, but is verified here by fitting the data. 
From |33]|, it follows that the rotation curves are self- 
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FIG. 4. Eight universal rotation curves with different ab- 
solute magnitudes Mi introduced in (33). The radial coor- 
dinate x in units of the optical radius x opt which encom- 
passes 83% of the total integrated light and the velocity v 
in units of v(x op t)- For the halo (long-dashed curve), the ini- 
tial values are: cr(0) = 3.8, 3.5, 3.3, 3.1, 2.7, 2.4, 1.9, 1.6 x 1(T 6 
of Mi = -18.29 for the first and Mi = -23.08 for the last 
value. 



similar, i.e. only one parameter, the luminosity or the lim- 
iting rotation velocity, establishes completely the proper- 
ties of the halo far from the luminous matter region. This 
property is also revealed in our model where asymptot- 
ically only the parameter A det erm ines the halo; cf. the 

Moreover, we can 



Newtonian solutions of Section III 



connect this one free parameter with the amplitude of 
the scalar field at the center. 

With the improvements of measuring rotation curves, 
HI rotation curve data have been found, investigating 
higher distances from the galaxy center. For an exponen- 
tially thin matter disk with a surface density T,(x/xq) = 
So exp(— x/xq), where xq is the disk scale length and So 
some constant at xq, the contribution to the rotation ve- 
locity can be determined by fill (cf. also H] ) 
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FIG. 5. Dwarf galaxies: Rotation curve fits for NGC1560, 
UGC2259, DD0154, and DDO170 with halo (long-dashed), 
stars (dashed), and HI gas (short-dashed). The velocity v is 
measured in [km/s] and the radial coordinate r in [kpc]. The 
parameters for these fits are summarized in Table 



(24) 



where /„ and K n are modified Bessel functions. Equa- 
tion ( p2| ) is actually the approximation of this exact re- 
sult. This formula shall be used in the following for the 
contribution of stars and gas in individual galaxies, with 
specific constants ro and So- The gas part contributes 
with a summand Vg as under the square root in (p3|). 

In some galaxies, a clear bulge can be read off the light 
curve. For these cases, we use the formula from Kent 
p5| in order to calculate the circular velocities from the 
observed surface density abuige (r) 



V b ulge( r ) = / C,a bu lge{C,)dC + 



K 

2nr 



arcsm — 



(VbulgeiOdC ■ (25) 



Corrections for flattening of bulges are ignored by this 
formula. 

The galaxies we choose belong to a sample of 1 1 galax- 
ies of different Hubble types and absolute magnitudes 
which fulfill several strong requirements p6| ]. All rota- 
tion curve data are measured in the 21-cm line of neutral 
hydrogen, so that the gas distribution extends far beyond 
the optical disc (at least 8 scale lengths) and the necessity 
of dark matter is becoming obvious. Isolation of galax- 
ies are another constraint so that perturbative effects of 
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FIG. 6. Dwarf galaxy: Rotation curve fit for NGC3109 
with halo (long-dashed), stars (dashed), and HI gas 
(short-dashed). The velocity v is measured in [km/s] and 
the radial coordinate r in [kpc]. The parameters for this fit 
are summarized in Table 0. 



nearby situated galaxies are negligible. Besides of that, 
high quality data are required. 

For the calculation of the neutral hydrogen gas for the 
galaxies of this sample, a further remark is in order. In 
[BjJ , it was shown that the HI data can be decomposed 
into a sum of two or three exponential disks. We shall 
use the results of this paper. We start now fitting dwarf 
galaxies. 

TABLE I. Data for rotation curve fits: luminosity L, disk 
scale length ro, scalar field frequency ui, mass-to-light ratio 
for disk and bulge, initial value of the scalar field a. In all 
cases, we have used v(0) = —5 x 10~ 6 . The parameters for 
the decomposed HI gas surface densities can be found in |j7| ; 
we have multiplied the total NGC3109 HI mass by a factor 
of 1.67 because of some missing 21 cm line flux at the VLA 
where the NGC3109 observations were made [HQ. [Hubble 
parameter Ho = 75 km/(s Mpc)]. 
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Dwarf galaxies have the characteristic of very low lumi- 
nosities p0| . Following empirical models, their rotation 
curves vary from the ones of spiral galaxies in such a way 
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FIG. 7. Spiral galaxies: Rotation curve fits for NGC2403, 
NGC2903, NGC3198, and NGC7331 with halo (long-dashed), 
stars (dashed), and HI gas (short-dashed); NGC7331 has also 
a bulge included which is truncated at 6.2kpc The ve- 

locity v is measured in [km/s] and the radial coordinate r in 
[kpc] . The parameters for this fit are summarized in Table [l| 



that instead of the Hernquist profile |3Jj] an isothermal 
density profile has to be applied |30| ; cf . also |3j| . From 
the rotation curves, it can be derived that the dark mat- 
ter density near the center is almost constant, i.e. dark 
matter has a core in these galaxies. Our Newtonian so- 
lutions for the massless complex scalar field reveal a con- 
stant density near the center, so it is not surprisingly 
that good fits for five dwarf galaxies from the Bcgeman et 
al. sample have been found (Fig. |B| and ||) 40 - 43|] . In all 
cases, it is recognizable that the dark matter halo dom- 
inates the luminous parts of the galaxies. Especially re- 
markable is the maximum of the DDOf 54 data which can 
be matched perfectly by our dark matter halo; hence, de- 
creasing rotation curve data, at the end of observational 
resolution, can be explained by using only a dominating 
dark matter component. The prediction of our model 
is to find oscillations in rotation curve data at high dis- 
tances from the center of galaxies. 

Our results for the sample of spiral galaxies are sum- 
marized in Fig. ^, cf. In all cases, we are able to 
produce good fits. Models using an approximate isother- 
mal sphere or a modified Newtonian dynamics (MOND) 
were also applied for these galaxies, e.g. [ (45pq . Confor- 
mal gravitation theory has been explored in |27| 47],|37],|^| ; 
the disadvantage of that model is its increasing rotation 
curves. 

By building universal rotation curves, a smoothing 
method is used [[33). Important informations about the 
distribution of the dark matter within individual rotation 
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FIG. 8. Spiral galaxy with oscillating rotation curves: 
NGC2998 with dominating halo (long-dashed), stars 
(short-dashed), HI gas (dotted), and bulge (dashed). The 
velocity v is given in [km/s] and the radial coordinate r in 
units of [kpc]. In order to find a good correspondence with 
the data it was necessary to have a dominating halo. The 
extrema of data and halo are at about the same values of r. 
The bulge is truncated at 2.9kpc; the data for the bulge are 
taken from [Q, but notice a different value for the Hubble 
constant there. The HI gas contribution |^(| was decomposed 
into (TBI = 55exp(-r/12kpc) - 49 exp(-r/9kpc) [M /pc 2 ]. 
For other parameters, see Table [j]. 



curves could be lost by this method. One can find an os- 
cillating behavior within the optical rotation curves of 
some galaxies, for example NGC2998 [E3]; cf. 'ripples' in 
the light curve profile [B5J . The explanation is that across 
a spiral arm, positive velocity gradients are observed; the 
velocity decrease from the outer edge of one arm to the 
inner edge of the next arm is in some galaxies faster than 
Keplerian what is taken as compelling evidence for non- 
circular velocities. In Fig |^, we fit the rotation curve 
data from Rubin et al. |49|] . We have taken into account 
four components: A truncated bulge, a thin disk of stars, 
an HI component, and a dominating dark halo. What 
one can recognize is that the data are not fit too well but 
that the maxima and minima of the fit and the data are 
at the same place. The dark halo has to be the main 
part in order to find a good fit. (The partly non-smooth 
behavior of the curves results from numerical problems 
because of the scaling of our solution by the scalar field 
frequency to.) Further investigations will show whether 
cylindrically or axisymmetric solutions can improve the 
fit or rule out our model. NGC2998 belongs to a sample 
of further 22 spiral galaxies for which the MOND model 
was applied |50|. Oscillations ('wiggles') have been found 
also in Hcv rotation curves [pit]. 
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V. SCALING OF THE SOLUTIONS AND THE 
MASS DISTRIBUTION 



(29) 



As we have seen from the redefined quantity x — wr, 
the frequency uj scales the physical dimension of the so- 
lution. Additionally, the mass is scaled with X/u> and 
the energy density with uj 2 . Table || shows for different 
values of oj the mass and the density at different radii. 

There exist two ways how one can interpret the values 
of Tabic [n]. The last row (uj = 10~ 21 /cm) shows a halo 
having a central density of 10 _21 g/cm 3 ; at x — 20 which 
means r — 20kLy the density decreases to a value of 
10 _23 g/cm 3 and within a sphere of this radius the halo 
has a mass of 10 43 g= 10 10 M Q . 

But one can use Table [l| also in the vertical direction. 
A solution with uj = 1 /cm (first row) has within a sphere 
of radius 20cm a mass of 10 24 g, within a sphere of 20 x 
10 15 cm (20 times the diameter of our solar system; second 
row) a mass of 10 37 g, within a sphere of 20Ly (third row) 
a mass of 10 40 g, and so on. The same procedure is valid 
for the density columns. The reason why one can do this 
is that one has a uj independence of the density and the 
mass at high radii, but a uj dependence for small radii. 
For the mass dl^), we find 



MM = —A 2 

K 



cos(2ujr) — 1 
2uj^ 



which goes over for small r values into 



M (r) oc —A 2 (-uj 2 r 3 --uj 4 r 5 
k I 3 45 



(26) 



(27) 



The same is valid for the density p. This result shows 
that the parameter A determines the behavior of the halo 
at high distances from the center (namely the limiting 
orbital velocity), but near the center also the frequency 
of the scalar field influences the behavior of the mass and 
the density. Rewriting Eq. (Era) into 



2uj 2 r 



nM(r) 

47T 



A 2 r 



cos(2wr) 



1 



(28) 



we immediately see that the left hand side has to be non- 
positive from which we can find the inequality 

TABLE II. Order-of-magnitude estimation for A = 10~ 3 . 
We put the frequency \juj = 1cm, 10 15 cm (about the radius 



of our solar system), 10 cm = lLy, and 



lkLy. 



Further: uj in units of eV; density p at the center and at a 
distance of x — 20; the mass within a sphere of radius x — 20. 



at an arbitrary radius R. This inequality states that at 
every radius R, the relation between the Schwarzschild 
radius and the radius of the object is smaller than the 
square of the central value of the scalar field which is 
10~ 6 for a galaxy with v max = 300km/s. 

Let us, for example, take the value uj = 10~ 21 /cm of 
the last row of Table H. To calculate the mass and the 
density for smaller radii, we have to take the formulae 
@ and @. The result can be seen in Table fill 
recognize in comparison with the values of Table 
stantial smaller values. 



The order of magnitude values of Table III for the mass 
and the density show also a possibility of a hierarchy of 
objects. Recently ]52| , it has been detected that the dark 
matter distribution in the Fornax cluster of galaxies is a 
mixture of two distinct components on different scales 
(not necessary distinct matter components). One part 
is associated with the cluster and the second smaller one 
with the galaxy NGC1399. This observation was done by 
looking into the X-ray-emitting diffuse plasma of tem- 
peratures 10 7 -10 8 K which is widely found in elliptical 
galaxies and clusters of galaxies. This gas is likely to be 
in hydrostatic equilibrium, where the thermal pressure is 
balanced by the gravity of the cluster. In this way, the X- 
ray emission traces the distributions of the dark matter 
in regions of higher density where more gas is trapped. 
This observation is the first evidence for a hierarchical 
nature of the dark matter distribution. Future investiga- 
tions will show whether two of our dark matter halos can 
be combined: one for the cluster, the other for the single 
galaxies. 



VI. ENERGY DENSITY, TANGENTIAL 
PRESSURE, AND REGULARITY 

An investigation of the energy density p (i.e. the radial 
pressure p r ) shows a terraced decrease (Fig. ^|). Fur- 
thermore, we find the relation | p \ > \ p± , where the 
tangential pressure p± oscillates sinusoidally around zero. 
From the equation of state, one recognizes that at each 
extremum of the scalar field a the difference between ra- 
dial and tangential pressure vanishes. Hence, at constant 
spheres, the equation of state p = p of a stiff fluid arises. 

TABLE III. Order-of-magnitude es- 

timation for uj = 10~ 21 /cm. In comparison with Table [q| 
we find smaller densities and masses at smaller distances. 



ui [1/cm] 


uj [eV] 


p(0) [g/cm 3 ] 


p(20) [g/cm 3 ] 


M(20) [g] 


distance to center 


P [g/cm 3 ] 


M [g] 


1 


10~ b 


10 2i 


10 iu 


10 24 


20 cm 


10~ 2i 


IO" 17 


io- 15 


10~ 20 


1Q -9 


io- 11 


10 37 


20 IO -3 Ly 


IO" 21 


10 28 


1Q -18 


1Q -23 


io- 15 


io- 17 


10 40 


20 Ly 


IO" 21 


10 37 


io- 21 


1Q -26 


io- 21 


1Q -23 


10 43 


20 kLy 


1Q -23 


10 43 
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TABLE IV. This table shows the radius R at which our 
dark matter halo reaches the density p. v is the limiting 
orbital velocity. 



300 



FIG. 9. The energy density p (solid curve) and the tan- 
gential pressure p± (dotted curve) [both in units of w 2 /k] for 
the interval of x — 200 up to x — 300. One recognizes very 
clearly the terraced decrease of the energy density. The mat- 
ter behaves like stiff matter at the saddle points of p. [Initial 
values: <r(0) = 1 and u(0) = 1] 



We check numerically the quadratic curvature invari- 
ant (Kretschmann scalar) Rf nv = R Xatlu R\ a ^u revealing 
its global regularity for our solutions |20| . The same was 
found for the two invariants R^R^ of the Ricci tensor 
and R? of the curvature scalar. The invariants of the 
irreducible decomposition of Riemann's curvature tensor 
show also no singularity [^0| . Hence our solutions possess 
no physically relevant singularity. 



VII. CONSTRUCTION OF A SURFACE 



<r(0) 


v [km/s] 


P fe/cm 3 ] 


R [Mpc] 


10~ 3 


300 


10" za g/cm ii 


2.3 






10" 34 g/cm 3 


830 


4 x 10~ 4 


120 


lO-^g/W 


0.9 






10- 34 g/cm 3 


330 



and we find w(2)/m = 0.3203. Within the surface, the 
metric potentials of the massive scalar field go over to 
the Schwarzschild metric. Hence, we have an object with 
a finite mass. In the next Section, we shall show another 
possibility to stop the halo. 



VIII. COVER FOR THE HALO 

The energy density p of the Newtonian solution de- 
creases in leading order as A 2 /x 2 . Therefore, we can 
come to the conclusion that the energy density will even- 
tually attain an hydrostatic equilibrium with another 
matter form at the same pressure value, e.g. the cos- 
mic microwave background radiation (CMBR), which is 
about 10 _34 g/cm 3 (seas of neutrinos and supersymmet- 
ric particles would increase this density). Alternatively, 
a cosmological constant A of the order of magnitude 
10 _29 g/cm 3 would produce much smaller halos; for es- 
timates of A see p3| . From Table IV , we derive a radius 



of about 10 6 pc for A and about 10 8 pc for the CMBR 
value (cf. the discussion in p4j). Recent observations 
show halo radii of about 400kpc pl| . Outside the halo, 
within the spacetime with a constant density, we have a 
Schwarzschild-de-Sitter metric. 



A massive complex scalar field with mass m can be 
used for constructing a surface of our halo. This massive 
scalar field has an exponentially decrease and produces a 
finite mass. We have two different coordinates x\ = mr\ 
for the surface with the massive complex scalar field and 
X2 = ujr2 in the interior with the massless complex scalar 
field. At the connection we require x\ = x%, hence 
Lu/m = n /r% ■ We need an equilibrium of the radial pres- 
sures and continuity of the metric potentials. Because 
of the different scalings in the two regions, we find for 
the radial pressures w(2)/m = \/p r (l) /jv(2) [argument 
(1), e.g. in w(l), describes the model of massive complex 
scalar field and argument (2) that of the massless one], 
i.e., if we have the equality of the metric potentials, the 
relation between the pressures gives us the relation be- 
tween the frequency of the massless scalar field and the 
mass of the massive scalar field. It was possible to find 
numerically such solutions, for example: Initial values 
for the interior region: ct(0) = 0.2333, i/(0) = -0.4055; 
for the surface region: a(0) = 0.1, i/(0) = -0.14035418; 
equality for the metric potentials was looked for at x = 3 



IX. GRAVITATIONAL REDSHIFT AND 
ROTATION VELOCITY 

By increasing the initial value of the scalar field <r(0), 
the Newtonian description breaks down, i.e. the formulae 
from Section III are no longer valid. The metric poten- 
tials v and A are not constant and the full nonlincarity 
of Einstein's theory has to be taken into account. 

The mass and the density function change their behav- 
ior smoothly. They still increase together with cr(0) 2 but 
one recognizes an additional decrease of the x dependent 
terms; Table |v| shows the deviation from Newtonian be- 
havior for the mass function with increasing <r(0). This 
means that one can still use the order of magnitudes of 
Table || (which was for <r(0) = 10 3 ) but one has to mul- 
tiply the values of the mass and the density by 10 4 for 
cr(0) = 0.1 and by 10 6 for <r(0) = 1, for example. The 
deviation of the formula is so eminent that one has to use 
numerical calculations. 
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TABLE V. Deviation for M/a 2 {0) from the Newtonian for- 
mula (^) with increasing <j(0); the mass still increases with 
<j(0). The values for x = 0.1 have to be multiplied by 10~ 3 . 



TABLE 
redshift function 2 



VI. 



The 



(a;) = exp\u(x) — v(Q)] — 1 for solutions 
with different initial values of the scalar field (f(0) = in 



X 


Newtonian 


<r(0) = 0.1 


0.5 


1.0 


1.5 


each case). 














<x(0) = 0.25 












0.1 


0.33288 


0.33287 


0.33255 


0.33156 


0.32991 


z 


n <{ 

u . o 




i 

i . 


1 9^ 


1 <{ 


1.0 


0.291 


0.290 


0.267 


0.212 


0.154 


8(10) 


0.16 


0.66 


1.58 


3.16 


5.88 


9.52 


2.0 


1.586 


1.571 


1.261 


0.747 


0.422 


z(50) 


0.25 


0.97 


2.13 


4.04 


7.83 


16.17 


3.0 


2.993 


2.955 


2.235 


1.193 


0.633 


2(100) 


0.29 


1.08 


2.31 


4.29 


8.29 


17.57 


4.0 


3.856 


3.798 


2.799 


1.513 


0.811 


z (200) 


0.32 


1.18 


2.47 


4.50 


8.66 


18.65 


5.0 


4.816 


4.715 


3.166 


1.739 


0.972 


2(300) 


0.35 


1.24 


2.56 


4.62 


8.86 


19.17 














z (400) 


0.36 


1.28 


2.61 


4.70 


8.98 


19.52 


























2(500) 


0.37 


1.31 


2.66 


4.76 


9.08 


19.77 



The strong curvature connected with these general rel- 
ativistic (GR) solutions produces also a gravitational red- 
shift z which, for a static spherically symmetric mass dis- 
tribution, is given by 



[v(A 1 )-v(A )]/2 



1 



(30) 



where Aq corresponds to the place of the emitter and 
A\ to that of the receiver, both in rest. Table VI shows 
results for different values of cr(0) up to x = 500. The 
redshift increases with growing a(Q) and with growing ra- 
dius. Redshift values reached by a(n emitting) gas cloud 
at the center of the halo can achieve very high values; we 
have not found numerically any limits. The z values for 
different x show the gravitational re dshi ft which appears 
if our GR halo ends at x. In Table VII, we assume that 



the gas clouds are in rest within the gravitational halo 
potential (more likely they circulate around the center 
so that an additional Doppler redshift has to be taken 
into account). Then, the redshifts can be interpreted 
as emission or absorption lines of the corresponding gas 
clouds. An interpretation of both Tables could be (for 
<t(0) = 1.0): Some radiation comes from the center and 
is gravitationally redshifted by z = 4.76, some part of 
the radiation is absorbed from a gas cloud at a distance 
x = 1 (z = 3.398), another part at x = 100 (z = 0.087), 
and finally the halo ends at x = 500 so that gravita- 
tional redshift does not take any further effect. Near 
the center (x = 1), we recognize a steep decrease of red- 
shift values (cf. Table VI and VII). The scale of such 
a halo is again defined by the frequency u>. If one as- 
sumes uj = 10 _18 /cm, then the halo has a diameter of 
r = 500Ly; hence, not very far away from the source, 
about r = 400Ly, one has already small gravitational 
redshift values. 

How fast does a mass rotate in such a spacetime? Fig- 
ure [To] shows a numerically calculated solution. We see 
that velocities of more than 10 5 km/s are reached. This 
means that 6% of the rest mass energy is stored in ki- 
netic energy. If we assume that per year a mass of 1M© 
transfers its kinetic energy into radiation then we find a 
luminosity of 10 44 erg/s. 



TABLE VII. Redshift 

values for absorption lines zix) = \J exp[^(500) — v(x)] — 1 
for different initial values of the scalar field (y (0) = in each 
case). We assume that the halo ends at x = 500. 



z 


o-(0) = 0.25 


0.5 


0.75 


1. 


2(1) 


0.350 


1.147 


2.124 


3.398 


2(5) 


0.229 


0.536 


0.630 


0.669 


2(10) 


0.183 


0.386 


0.417 


0.383 


2(50) 


0.096 


0.170 


0.168 


0.141 


2(100) 


0.064 


0.108 


0.105 


0.087 


2(200) 


0.034 


0.056 


0.054 


0.045 


2(300) 


0.019 


0.030 


0.029 


0.024 


2(400) 


0.008 


0.012 


0.012 


0.010 


2(500) 


0. 


0. 


0. 


0. 



X. DISCUSSION 

We have shown that we were able to fit the universal 
rotation curves of Persic and Salucci and of a selection 
of dwarf and spiral galaxies which are especially suitable 
to find out characteristics of isolated dark matter halos. 
In one case, we could match oscillations of rotation curve 
data and our model. 

We have to use Newtonian solutions and to fix two 
parameters. There is the central amplitude of the scalar 
field, the parameter A, which determines the limiting 
orbital velocity. The second parameter, the frequency 
of the scalar field, uj, varies the values of the mass and 
the density near the center of the galaxies. 

Because our solution consists of a condensed 'star'-like 
object (with increasing mass) also the energy contribu- 
tion of the radial pressure p r plays a role for the rotation 
velocity. We found that this part contributes the same or- 
der of magnitude as the Newtonian part coming from the 
'normal' energy density p. Therefore, general relativity is 
just the correct theory for the dark matter part provided 
a compact object with internal pressure is present. 

Another issue is the physical basis of the Tully-Fisher 
relation. It is obvious that in a deeper gravitational po- 
tential, i.e. more dark matter, also more luminous matter 
can find place. Our best fits for the rotation curves are 
just found if we add for low luminosity galaxies an al- 







1.0 , 1.7 




5 10 15 SO 5 10 15 20 

X X 

FIG. 10. The scalar field a, the metric potentials e" and 
e A , and the rotation velocity v (in units of [km/s]) for high 
relativistic initial values of <r(0) = 1 and v(0) = —3.5. Near 
the center, particles in circular orbits have velocities of more 
than one third of the velocity of light. 



tic boson stars were used as dark halos. In the latter 
case, the contribution of the radial pressure to the ro- 
tation velocity was wrongly neglected. It is well-known 
that excited boson stars are unstable against small ra- 
dial perturbations (see e.g. p3|); for a recent numerical 
investigation of instability of excited boson stars, see |i0) . 
Under the assumption that boson stars are transparent 
(as the solutions of this paper), maximal gravitational 
redshift values of about 0.7 have been revealed pn{ ]. In 
p2| , it was shown numerically and analytically that the 
Newtonian solutions of our model are stable. 
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most dominating scalar field halo and vice versa. The 
contribution of our dark matter form is determined by 
the parameter A the values of which lie in the order of 
magnitudes 10 -4 and 10~ 3 . In a future paper, we will 
investigated why such values are preferably found. 

In a recent paper [|55| , the reason for the Tully-Fisher 
relation is divided into two parts: (a) a mass-to-light ra- 
tio of the luminous matter and (b) a relation between the 
luminous mass and the rotation velocity. It is shown that 
there is a smaller scatter in part (b) than (a), because (a) 
depends on the present star formation rate. The conclu- 
sion is that part (b) in combination with a well-behaved 
relation between luminous and dark matter (producing 
flat and smooth rotation curves, i.e. the 'conspiracy') is 
the physical basis of the Tully-Fisher relation. The con- 
spiracy maintains a luminous-mass-rotation-velocity re- 
lation with a slope of 4. The issue is whether one can 
reveal a reason from our model. We will investigate this 
point in a future work. 

Providing this halo model is realized in nature and one 
finds a radius of the dark matter halos of galaxies, then 
this radius could be a hint for a cosmological constant 
and its order of magnitude. 

During the completion of this paper investigations with 
massive scalar particles were carried out. In [56 57|, an 



isothermal ideal Bose gas which is degenerated in the cen- 
ter of galaxies was used to fit rotation curves of 36 spiral 
galaxies. The authors found the best fits with a boson 
particle mass of about 60eV. In |^8| excited Newtonian 
boson star solutions and in |p9(| excited general relativis- 
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